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Abstract 

We calculate the next-to-leading order(NLO) quantum chromodynamics(QCD) corrections to 
the inclusive process of Zq —>• B* + c + b under the non-relativistic QCD(NRQCD) factorization 
scheme. Technical details about contributions from vector and axial-vector currents in dimensional 
regularization scheme are discussed. Numerical calculation shows that the NLO correction enhances 
the leading-order decay width by about 50%, and the dependence on renormalization scale g is 
reduced. The uncertainties induced by quark masses and the renormalization scale g are also 
analyzed. 
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I. INTRODUCTION 


Due to its unique nature in the family of mesons, B c system attracts wide attention of ex¬ 
periment and theory. The study of B c mesons may deepen our understanding of the Standard 
Model(SM) and the effective theory, non-relativistic quantum chromodynamics(NRQCD) 


1]. B c meson was first discovered at the Fermilab TEVATRON 0, and its excited state 


B C {2S ) [3| was recently observed by ATLAS Collaboration. 


S 


For B c meson direct production, various investigations had been carried out. Refs 


studied in detail the B c hardroproduction, in which the P-wave [5j] and color-octet 


■ft 


6 ] 


contributions, as well as the “intrinsic heavy quark mechanism” 171 were taken into account. 

n q 

Refs. [8|, [9[ discussed the B r production through the fragmentation scheme, and a Monte- 
Carlo simulation program 10] for the B c hadroproduction was also produced. 


Apart from the direct production, B c meson indirect production is also interesting, which 
may inform us not only the nature of B c meson, but also the characters of its parent particles. 
B c production through the top quark decays was discussed in Refs. (ill . 0 , through W ± 


decays was calculated in Ref. 


In the work of Ref. 
calculated. 


3 . 


n 


13|, and through Z° decays was analyzed in Refs. 0flll5l. 
the NLO QCD corrections to Z° —$■ R c ( 1 S'o) + c + b process was 


At the LHC and International Linear Collider(ILC) 17f, or other forms of Z factories, 
the Z° boson is and will be copiously produced. The Z° production cross section at the 
LHC is about 34nb 18|, and at the ILC, e.g., the cross section will be about 30nb while 
collider runs at the Z° pole energy 0], Given the colliders’ luminosity to be 10 34 cnU 2 s _1 ~ 
10 8 nb _1 /year, there will be ~ 10 9 Z° events being produced per year at the LHC and ILC. 
Therefore, to study the B c production in Z° decays is worthwhile and meaningful. And, 
for this aim, since the B* will almost completely decay to scalar B c and the NLO QCD 
correction in heavy quarkonium energy region is large, in some cases even huge, in this work 
we calculate the NLO QCD corrections for Z° —> B^^Si) + c + b process. 


The rest of the paper is organized as follows. In section II we recalculate the Z° to 
B* process at the Born level. In section III, the NLO virtual and real QCD corrections 
to the leading order result are evaluated. Section IV presents some technical details of 
the calculation. In section V, the numerical evaluation for the concerned decay process is 






















FIG. 1: The leading order Feynman diagrams for B* production in Z° decays, 
performed at NLO accuracy. The last section is for summary and conclusions. 


II. LEADING ORDER CONTRIBUTION 


The Z°(k) -A B*(p 0 ) + c(ps) + b(pe) process starts from a 2 order, the Born level, and the 
corresponding four independent Feynman diagrams are shown in FigJTl The decay width 
may be expressed in a standard form, i.e. 

dT Born = V |-M Born 1 2 dPS 3 . (1) 

2 mz 3 

Here, Y me a n s summing over the polarizations and colors of final particles, 1/3 comes from 
the spin average of initial state, and dPS 3 stands for the three-body phase space of final 
states, which can be explicitly expressed as 


dPS 3 


1 

2 7 7r 3 ?n| 


ds 2 dsi , 


( 2 ) 


where = (p 0 + P 5) 2 = [k - p 6 ) 2 and s 2 = (ps + Pe) 2 : 
bounds for and s 2 can be found in our previous work 
particle physics. 


(k — p 0 ) 2 . The upper and lower 


16], or any standard text book of 


The amplitude M Born = M a +M b +M c +M d can be readily obtained from Fig ]T] according 
to Feynman rules: 


M a 


c x u(p 6 )y* 


v(P4)u(p 3 ) 

(P4 + Pe) 2 


min 


7 5 ) 


(- ^4~ ^5~ A) + ™c 

(-P4 - Ps - Pe) 2 - m 2 


7/XPs) , 


M b 


C x u(p 6 ) 7 m 


(j4+ #5+ A) + m b 

(P3 +P5 + Pe) 2 - m 2 b 


Ak){ 7 5 


n) 


n(p 4 )n(p 3 ) 
(P3 +P5) 2 


7m^(P5) 


J 


-M c 


C x n(p 6 )7 m 


u(p4)n(p 3 ) 
(P4 + Pe) 2 


7m 


(^ 3 + ^4+ j4) + 

(P3 +P4 + Pe) 2 - "i 2 


/(fc)(T c - 7 >(p 5 ) 


J 


3 











M d 


C x u(p 6 )ftk)( 7 5 


n) 


(- fa- fa- lf 5 ) + m b v(p 4 )u(p 3 ) 
{-P3 ~Pa -ph) 2 -ml' (P3+P 5 ) 2 ' fl 


( 3 ) 


Here, the constant C = with 6\y being the Weinberg angle, e(k) is the polarization 

vector of Z° boson, T c = (1 — | sin 2 9w) and T b — (1 — | sin 2 d\y). For c and b constituent 
quarks hadronization to B* meson, the following projection operator is employed 0 


v(p 4 )u(p 3 ) -a 


^3 5l (0) 

2 V m ^ 


f(Po)(tfo + m B *) 



( 4 ) 


where e*(p 0 ) is the polarization vector of B* meson with p 0 = p 3 + p 4 , m B * = m c + m b , 
l c stands for the unit color matrix, and v I / 3s 1 (0), a nonperturbative parameter, is the 
Schrodinger wave function at the origin of the B* meson. 


III. NEXT-TO-LEADING ORDER CORRECTIONS 

The NLO QCD corrections to the Z°(k) —* B*(p 0 ) + c(p 5 ) +b(p e ) process contain virtual 
and real corrections, i.e. Y Virtua i and Y Rea i respectively, which are both at the order of a^. 
Typical Feynman diagrams which attribute to the virtual correction are presented in Figsj2j- 
[5] while those for the real correction are shown in FigjGj Note that in Figsj3]- |6]only those 
diagrams with Z —y cc vertex have been displayed, the remaining half can easily be obtained 
by exchanging the c and b quark lines. 


A. The Virtual Correction 


With virtual correction, the decay width can be formulated as 

1 1 


dr 


Virtual 


2 mz 3 


J2 2 MM V irtualM* Bo rn)dYS 3 . 


( 5 ) 


In Re(M V irtualM* Born ), both ultraviolet(UV) and infrared(IR) divergences exist. The con¬ 
ventional dimensional regularization scheme wiht D = 4 — 2e is adopted to regularize them. 
There are also Coulomb divergences, which in this work are factorized out through the 
threshold expansion technique 21] and then attributed to the bound state wave function. 
In the calculation, the NRQCD short distance coefficients are obtained by matching to the 
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FIG. 2: The typical counter-term Feynman diagrams corresponding to FigJTIa). 


full QCD result stemming solely from the hard interaction region, and other regions give no 
contribution. 

According to the power counting rule, the UV divergences exist merely in self-energy 
and triangle diagrams, which are canceled by counter terms(CT). The Figj2] contains 5 
typical CT diagrams corresponding to the FigUJa), and hence there should be other 15 CT 
diagrams not shown. Of the 20 CTs, the renormalization constants include Z 2 , Z 3l Z m , and 
Z g , corresponding to the renormalizations of quark field, gluon field, quark mass, and strong 
coupling constant a s , respectively. In practice, the terms related to Z 3 vanish, e.g. those 
contribute to Z 3 in Figj2j, the ( 03 ), ( 04 ) and ( 05 ), cancel with each other. In the calculation, 
Z g is defined in the modified-minimal-subtraction(MS) scheme, while for Z 2 and Z m we take 
the on-shell(OS). I 11 the end of the day we have 


au 


5Z^ b = -C F — 
An 

a < 


1 2 

+ —-3 ln(m 2 ) + 4 


L e uv 


-ir 


5Z™ = 


-3 C 


An 


L e uv 


ln(m 2 ) + 


+ Otf) , 
+ , 


f 7 MS _ PO a s 

9 2 An 


L e UV 


- ln(/i 2 ) 


+ O(o2) 


( 6 ) 


Here, 1 / t UV ( JR) = 1/e — 7 ^ + ln( 47 r / u 2 ), /i is the renormalization scale, the mass m in 
5Z° S and 5Z^ S stands for m c and accordingly, /3 0 = {11/Z)Ca — (4/3 )TfUf is the 
one-loop coefficient of the QCD beta function with rif = 5, the number of active quarks 
in our calculation. Eventually, the UV divergences appearing in CTs eliminate all UV 
divergences in self-energy and triangle diagrams, leading to a UV-free result. The remaining 
IR divergences in CTs have the following form 

2 Cpct 


dr & = 


-dr 


7 re 


Born • 


(7) 


IR 


In one-loop Feynman diagrams, the IR divergences involve triangle, box and pentagon 
diagrams. Of the triangle diagrams in FigQ] only two have IR divergences, i.e. TriangleNg 
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FIG. 3: Half of the self-energy diagrams in the virtual correction. 



FIG. 4: Half of the triangle diagrams in the virtual correction. 

and TriangleN 10 . Of the diagrams in Figj5j BoxN 7 has no IR divergences, BoxN ] and 
PentagonNx have both Coulomb singularities and IR divergences, and the remaining other 
diagrams have only the IR divergences. We find that the cancelation of IR divergences in 
Figs. (Bl and (J5]) goes as follows: 

• Combinations of TriangleNg + B 0 XN 2 , BoxNs + BoxNg are IR finite; 

• Combination of TriangleN 10 + BoxN 3 + PentagonN 2 are IR finite; 

• The remaining IR divergences lie in BoxN 1; BoxN 4 and PentagoiiN |: 

(I 













FIG. 5: Half of the box and pentagon diagrams in the virtual correction. 

• Diagrams with Z —» bb vertex have the same cancelation pattern as in above. 


B. The Real Correction 


For the real correction of the concerned process, in which the IR divergences exist, 
there are 24 Feynman diagrams and half of them are shown in FigjGl To regularize the 
IR singularities, the “Two cutoff phase space slicing method” has been employed 22], For 
diagrams contain IR divergences, we enforce a cut 6 on the energy of the radiational gluon, 
the pj. The gluon with energy < 5 is considered to be soft, while > <5 is treated as 
hard. Then, the decay width can be written as: 


dr Real — dr Rea i \ IR +dr Reai I iR-f r ee i 

dr M |„= dr™7”'‘ | p?<i +drgy“" i | p?>5 


(8) 


In FigU we find that 
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FIG. 6: Half of the Feynman diagrams in the real correction. 

• The first 4 diagrams, the RealNi^^, are IR free; 

• The combination of RealN 5 + RealN 6 , RealNy + RealNg has no IR singularity; 

• The remaining 4 diagrams, i.e. RealN 9jl0illjl2 , have IR divergences; 

• The diagrams with Z —» bb vertex behave the same as in above. 


Then the decay width in the soft sector dT^ aZ so -^ \ p o <s can be written as 


1 1 


drftr'' U<5= x dPs„ u. 


2 mz 3 

In the Eikonal approximation, the squared amplitude |A4 


I-Kea7 0/ T = -(C F 47ra s )\M B orn\ 2 X 


pi 


- 2 - 


IR—soft 1 2 j 

Real I reads 

P5-P6 


+ 


pi 


( 9 ) 


( 10 ) 


(p 5 • p 7 ) 2 (p 5 • p 7 )ip 6 ■ p 7 ) (p 6 • Pi) 2 . 

which contains all the IR singularities in the real correction, i.e. RealNg^opi,^ and the other 
4 diagrams with Z° — y bb vertex. The 4-body phase space in soft sector possesses the form 

cl 3 p 7 


dPS 


4 |soft 


= dPS, 


(27t) 3 2 p® 


l Py<6 


( 11 ) 
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Finally in the small S limit, the dr^ a / s °' P | p o <( 5 can be expressed as 


dr j£ai SO/ * U«5= dr Born — ^ ^-hi(<5 2 )^ ^1-+ finite terms ( 12 ) 

with x s = 1 — (^~ _ 4 mfe 2 m ^ c 2 ) T- Here, those 1/e^ involved terms in (TT2T) will cancel the IR 

singularities in CTs, i.e. (|7|) , and those remaining ones in the one-loop Feynman diagrams. 
Note that the ln(5 2 ) involved terms will be canceled by the 5-dependent terms in the hard 
sector. In the end, the final result will be IR finite. 


In the case of hard gluons, the decay width reads 


jt^IR— hard | 

Ui Real Ip7°>5- 


1 1 


2 mz 3 


s2> 


IR—hard |2 
Real I 


4 |hard 


(13) 


Here in |A4^ a / iarrf | 2 , the radiation gluon is considered to be hard, and the phase space 


dPS 4 | hard, can be expressed as 

2 \J(sy — m 2 c — ml) 2 — 4m;?ra 2 




i'2'K 


dPS 


4 |hard 


(47r) 


dp 0 ° J d cos 6 c 


d 0 c 


rP7- 


<-y+ 


X 


dp 7 ° 


d y + 


CP7+ 


n y+ 


'y- 


'P7- 


dp 7 ° 




d y 


(14) 


with 


P7° + = 


P V = 


s - 2y/sp 0 


2(y/s-p 0 ° - |po|) 


s - 2y / ipo° 

2(\/i — Po° + |po|) ’ 


y + = ~[(Vs - Po° - P7°) 2 - \Po\ 2 - (P7°) 2 + 2|po|p7°] 

y- = ~[(Vs - Po° - P7°) 2 ~\Po\ 2 - (P7°) 2 - 2|po|P7°] 
s 


(15) 


Here, y is a dimensionless parameter defined as y = (k—po—p 7 ) 2 /s, |po| = V ( Po°) 2 — 
and y/s = mz- 


The IR-free decay width with real correction can be formulated as 

dP-ReaZ | IR-free= ^^3 ^2 I^Rea/^T dPS 4 1 ( 16 ) 

where |A4j^~/ ree | 2 is the amplitudes squared without IR singularities, and dPS 4 has exactly 
the same form of dPS 4 |hard while 5 = 0. Finally, the sum of the soft and hard sectors, i.e. 
(fl2|) and ()T3j) . together with the IR-free part (TT6|) give the full real correction. 
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In the end, with real and virtual corrections, one can readily obtain the total decay width 
for the inclusive process Z°(k) —> B*(p 0 ) + c(p 5 ) + b(p 6 ) at the NLO accuracy of QCD, 

r total r Born “I - r Virtual T ^ Real T 0(cT s ) . (17) 


In (TT7T) . the decay width T tota i is UV and IR finite, and technical cut 5 independent as 
expected. 


IV. SOME TECHNICAL DETAILS IN THE CALCULATION 


In the conventional dimensional regularization scheme, the y 5 problem needs to be han¬ 
dled carefully, especially in the process which contains the axial-vector current. In this work, 
we adopt the scheme given in Ref. 23], where the following rules are followed. 

I. The anticommutation relations, i.e. {y'by 5 } = 0 and {y^y 1 '} = 2 g^ u . 


II. The cyclicity is forbidden in traces involving odd number of y 5 . When several diagrams 
contribute to one process, one should write down the amplitudes starting from the same 
vertex, named reading point. 


III. As a special case of rule II, in the anomalous axial-vector current situation, the 
reading point must be the axial-vector vertex in order to guarantee the conservation of the 
vector current. 


When applying these rules to our process Z°(k) —* B*(p 0 )+c(p 5 )+b(p e ), some conclusions 
are obtained. That is, in the virtual correction, the amplitudes squared have the following 
two structures, schematically shown in Fig® 

Structure 1: The two Z° — qq vertexes lie in one fermion trace. 

Structure 2: There are two fermion traces, each with a Z° — qq vertex and one involving 
triangle anomalous diagram, e.g. TriangleNn and TriangleNi 2 in Fig® 


For those whose amplitude squared belongs to ” Structure 1”, since the amplitude can be 
separated into vector part and the axial-vector part, amplitudes squared can be written as 

TVa ce\M Lmp M' Tr J = Tr[(MT^ + 

_ 'TVJ A Avector \ a vector 11 . rpr \ A axial—vector \ a axial—vector^] 

~ J- r 1 /W Loop JVl Tree \ + lI [ JVl Loop JVl Tree \ 

i r Pv.r \ Avector k a axial—vector^ , rpr a a axial—vector a a vector ti /i o\ 

+ 11 Loop M-Tree \ + 11 [M Loop ^ Tree J * 



Loop 


Tree 


non — Triangle 


Structure 1 


Structure 2 


FIG. 7: Typical structures of the amplitudes squared of virtual correction. 


In the last equation, though the first trace has no y 5 and the second trace has two 7 5 s, we can 
move them together by anticommutation relation and contract them to unit 1; for the other 
two terms each has one y 5 , and they will be canceled by their own complex conjugate terms. 
As a result, there is no need to evaluate the last two traces in (fT8|h and the calculation is 
hence greatly simplified. 

For the triangle anomalous diagrams case, or amplitudes squared in ” Structure 2”, as 
shown in Figj7]the amplitude squared reads 


AT Anomalous AT £j orn Trace [Ad Triangle\ TraC non—Triangle] 

_ r T r r A A vector . \ a axial—vector ] mr A A vector _i_ a a axial—vector ] 

II [y l Triangle * ^ ^Triangle J W *■ non—Triangle ' ^ ^non—Triangle! 

= T'AM^ie) Tr [M v :^ Triangle \ 

I rr r k / vector | m r k /axial—vector ] 

T r 'Triangle! W ' non—Triangle! 

, T r \ /axial—vector] m ,[ / /vector ] 

J.I [y 'Triangle J W 1 non—Triangle! 

i tv-T A /axial-vector 1 rp \ \ /axial—vector i /in\ 

+ ^[J^Triangle \ ^[-^non-Triangle! ' 


Here in the last equation, the first term will be canceled by the one whose triangle fermion 
loop is reversed, i.e. the vector currents of TriangleNn and TriangleNi 2 cancel with each 
other; the second and third terms will be canceled by their complex conjugate terms; the 
last term survives, and it is the only one we need to trace and contains one y 5 . In the final 
numerical calculation, the contribution coming from triangle anomalous diagrams turns out 
to be numerically insignificant. 

In the case of real correction, since all the amplitudes squared have the same structure 
as ” Structure 1”, one only needs to handle the Dirac traces without y 5 or with two y 5 s. 





In our calculation, the Mathematica package FeynArts 24] is used to generate the Feyn¬ 


man diagrams, FeynCalc 


25j and FeynCalcFormLink 26] are used to handle the algebraic 


trace manipulation, $ Apart 27], FIRE 28] together with codes written by ourselves are 


employed to reduce all the one-loop integrals into master-integrals(A 0 , B 0 , C 0 , D 0 ), and the 


LoopTools 


29] is employed to calculate the master-integrals numerically. The numerical 




integrations of the 3- and 4-body phase spaces are performed by VEGAS [30J. 


V. NUMERICAL RESULTS 


For the numerical calculation, the following input parameters are used 31]: 


m c = 1.5 ± 0.1 GeV , m;, = 4.9 ± 0.2 GeV , mz = 91.1876 GeV , rriw = 80.399 GeV , 


Gf 


sin 2 9 w = 0.2312 , g = 2\l2m w J-£= = 0.6531 , tf B .(0) = 


m 


= 0.3615 GeV 2 . (20) 


V2 ' ’ x/uF 

Here, Gf is the Fermi constant in weak interaction and -R(O) is the B* meson’s radial wave 
function at the origin, which value is estimated via potential model 32], The two-loop 
strong coupling of 

<**{») _ 1_ Pi In L 

4vr /3 0 L /3$L 2 

is employed in the calculation. In which L = /Aq CD ) with A qcd to be 214 MeV, and 

Pi = (34/3 )C\ — ACpTpnf — (20/3)C'aFf'^/ is the two-loop coefficient of the QCD beta 
function. 


In FigJH] the decay width r(/r) and the ratio T(/i) /T(2rrib) versus renormalization scale 
g are presented. The NLO corrections enhance the LO contribution to Z° —)■ B* + c + b 
when g > 4.6 GeV. We can see from Figj8]that the renormalization scale dependence of the 
decay width is reduced evidently. 

We adopt three typical renormalization scales in evaluation, and the corresponding values 
of running coupling constant are a s (/i = 2 mb) = 0.1768, a s (fi = 28 GeV) = 0.1423 and 
a s (/i = mz/ 2) = 0.1306, the decay widths are as presented in TabJH where the first errors 
are induced by m c = 1.5 ±0.1 GeV and the second ones are induced by mi = 4.9 ±0.2 GeV. 
Here, the K-faetor is defined as r r // o ° . In the calculation, when estimate the uncertainties 
induce by varying the m c , the mb is fixed and vice verse. 
















rW(keV) 


r(/i)/r(2m ft ) 




FIG. 8: The decay width T(^) and the ratio r(/r)/T(2m,&) versus renormalization scale for [i 
running from 2m c to mz/ 2. 

TABLE I: Decay widths T(/z)(keV) of Z° —>• B* +c + b at leading order and next-to-leading order. 
The first error are caused by m c = 1.5 ± 0.1 GeV and the second ones by mj = 4.9 ± 0.2 GeV at 
fj, = 2 mb, 28 GeV, mz! 2. 


r( M )(keV) 

LO 

NLO 

K-factor 

fj, = 2 mi, 

H = 28 

M = m z /2 

on cc —17.10+0.71 
0».0U +22 85 _o 71 

ro rjo —11.08+0.46 
oo.uz +14 80 _ 0 46 

a o 07—9.33+0.39 
^±0.0 1 -|_ 1 2.47—0.39 

118 77—23.11+0.37 
110 .1 ‘ +31.53+0.07 

n A 4n —18.25+0.45 
y4.4U_ | _ 24 .73_ 0 .22 

Q/i 16.33+0.44 

ou.uu+22 q 8 _o.26 

1 99 — 0.01 — 0.01 
1 -°°+0.01+0.01 

1 — 0.00 — 0.01 

1 - oo +0.01+0.01 

-1 70—0.00—0.00 
1 -*°+0.01+0.01 


The results of Table [D indicate that 

• The correction of NLO is significant, and the K-factor grows with the renormalization 
scale /i increase, yet it grows slower at high scale region; 

• When renormlization scale /j increases, the decay widths for both LO and NLO de¬ 
crease, yet both decrease slower at high scale region; 

• The uncertainties caused by varying m c are much larger than those induced by m&. 


To show the above results more clearly, in Figj9] we exhibit the LO(left) and NLO(right) 
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NLO 




FIG. 9: The decay widths of LO(left) and NLO(right) with different quark masses versus running 
renormalization scale //, from 2 m c to m z / 2. 

decay widths with different quark masses versus running renormalization scale (i. It is 
obvious that the decay widths are much more sensitive to m c then mb, the three lines of 
mb = 4.9(±0.2) are very close to each other at both LO and NLO. 

According to Ref. H and above calculation, the B c and b[*' production rates in Z° 
decays are readily obtained in LO and NLO, at scale /j, = 2 mb for illustration. That is 

Br LO (Z° ^B c + cb ) = Vl °W = 2 Q x 10 -5 

I z 

Br NLO (Z° ^B c + cb) = = 3.1 x 10" 5 , 

I z 

Br LO (Z° B* + cb) = Tl °( B ^ = 3.6 x 10~ 5 , 

l z 

Br NLO (Z° -)■ B* + cb) = Tnl o( B c) = 4_ 8 x 10 -s ? (22) 

I z 

where T z = 2.5 GeV is the total decay width of Z° boson. Because the spin-triplet state B* 
decays to the ground state B c with almost 100% rate, we obtain the total production rates 
for B c production in Z° decay, 

Btlo(Z 0 —» B c + X ) = 6.5 x 10- 5 , 

Br NLO {Z° ->■ B c + X) = 7.9 x 10" 5 . (23) 

Given about 10 9 or more Z° events being produced per year in future colliders, there would 
be sizable B c being produced. And the effects of the NLO corrections might be captured in 
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experimental observation. It is worth noting that the LO contributions from P-wave/color- 
octet are an order of magnitude smaller than the S-wave |li|, which only enhance the total 
production rates slightly. 


VI. SUMMARY AND CONCLUSIONS 

In this work we calculated the decay width of Z° —> B* + c + b inclusive process at 
the NLO accuracy in the framework of NRQCD. The calculation procedures for both LO 
and NLO corrections were presented. Our analyses for the vector and axial-vector currents 
contributions were performed in dimensional regularization scheme. The decay width and its 
uncertainties caused by varying quark masses, as well as the dependence on renormalization 
scale p were presented numerically. Supposing that there will be copious Z° produced in 
future colliders, our calculation together with previous work Q would be helpful for the 
precise study of B c physics, and might also tell how well the perturbative calculation and 
non-relativistic quark model work for B c system. 

According to the calculation, the NLO QCD correction to the inclusive process Z° —)■ 
B* + c + b is significant. We found that the renormalization scale p dependence of the decay 
width is depressed while the NLO correction is taken into account. When the scale p runs 
from 2 mb to mz/ 2, both values of T lo and Tnlo decrease, yet the ratio Ynlo/Ylo grows 
from 1.33 to 1.73, and the increasement trend slows down at high /i region. Moreover, the 
input parameter of quark mass m c has a quite large influence on the decay width, for both 
T lo and T nlO) at both low and high p regions as well. In contrast to m c , the uncertainties 
induced by rrii, are negligible. 
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